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We calculate the symmetry energy in cold dense matter both in the normal quark phase and 
in the 2-color superconductor (2SC) phase. For the normal phase, the thermodynamic potential 
is calculated by using hard dense loop (HDL) resummation to leading order, where the dominant 
contribution comes from the longitudinal gluon rest mass. The effect of gluonic interaction to the 
symmetry energy, obtained from the thermodynamic potential, was found to be small. In the 2SC 
phase, the non-perturbative BCS paring gives enhanced symmetry energy as the gapped states are 
forced to be in the common Fermi sea reducing the number of available quarks that can contribute 
to the asymmetry. We used high density effective field theory to estimate the contribution of gluon 
interaction to the symmetry energy. Among the gluon rest masses in 2SC phase, only the Meissner 
mass has iso-spin dependence although the magnitude is much smaller than the Debye mass. As 
the iso-spin dependence of gluon rest masses is even smaller than the case in the normal phase, we 
expect that the contribution of gluonic interaction to the symmetry energy in the 2SC phase will 
be minimal. The different value of symmetry energy in each phase will lead to different prediction 
for the particle yields in heavy ion collision experiment. 
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I. INTRODUCTION 


There is a world wide interest in symmetry energy as it provides a key to understanding the iso-spin asymmetric 
property of the nuclear matter. Furthermore, related topics range from rare isotopes at low energy regime to neutron 
star core at high energy regime ill, H]. The value of the symmetry energy at normal nuclear matter can be inferred 
from either the binding energy of semi empirical mass formula Q or experimental information from isotopes or heavy 
ion collision [2]. However, the value at higher density is still under debate. Some model calculations including iso-spin 
dependent interaction channel (NLpd) [j] predict a stiff symmetry energy as density becomes higher, which can cause 
less favourable condition for the neutron and subsequently decreased iso-spin density. 


The isospin asymmetric matter could persist longer in a heavy ion collision if the high density phase allows for a 
quark and hadron mixed phase. In such cases, it was found that the iso-spin density can remain high even if the 
hadronic part of the symmetry energy rises stiffly at high density 0-0- This will lead to the enhanced iso-spin rich 
resonances and subsequently to the increased iso-spin asymmetric decays. This means that the particle yields can be 
modified not only by the nuclear symmetry energy but also by the quark part. As representative studies, the MIT 
bag model was used to calculate the symmetry energy when the quark phase is normal 0, 6|, and additional color 
superconducting effect was taken into account through an NJL type model calculation [?']. Also the confined iso-spin 
density dependent mass (CIDDM) model was used in Ref. 0 to obtain the quark matter symmetry energy and related 
influence in the quark star. 


In the extremely dense condition, one should consider the equation of state for quark matter. The equation of 
state for non-interacting quarks with up to three flavors was first introduced in Ref. 9j in relation to the hypothetical 
quark stars. Due to asymptotic freedom, methods based on perturbative quantum chromodynamics (QCD) can be 
applied to calculate physical observables in extreme condition. If there exists some hard excitation scale in the matter, 
soft excitations can have infinite number of hard loop corrections that are not suppressed by the order of coupling 
and that can be resumed to an equivalent order 10-1121]. Considering the cold normal dense matter first, in-medium 
quark excitation scale is hard due to its large chemical potential, whereas the gluon excitation scale is mainly soft. 
Therefore, all equivalent gluonic excitations which contains hard dense quark loop corrections should be resumed [HI . 
Secondly, when temperature becomes even lower, all the fermions will be confined below their Fermi sea without 
large fluctuation. As the momentum fluctuation scales to a small value, most of the operators can be absorbed into 
the ener gy term or neglected. On the other hand, the four-fermion interaction with opposite momenta becomes 
relevant [14-17]. If the interaction is attractive, it is natural for the fermions to form a condensate which leads to 
superconductivity: BCS paring (10 . For QCD matter, as the quarks in the color anti-triplet channel are attractive, 
color superconductivity arises naturally with large gap size [19|-|23|. At the density regime where the quark-hadron 
mixed phase may exist, 2-color superconductivity (2SC) will be favored as there is a mismatch of Fermi momentum 
between light quarks and strange quark [20. 


In this work, we calculated the symmetry energy in cold dense matter using QCD. The density region being con¬ 
sidered lies between 3 to 5 times the nuclear matter density po = 0.16 fm~ 3 , which could be reached in heavy ion 
collision experiment. For the normal phase, we use the hard dense loop (HDL) resummation for the gluonic interac¬ 
tion. Thermodynamic potential is obtained to leading order in HDL perturbation, which is an extension/modifying of 
previously reported works [25l ~ l27j . The symmetry energy has been calculated from the thermodynamic potential. For 
color superconductor phase, 2SC phase is considered. Without perturbative gluonic interaction, we adopted thermo¬ 
dynamic potential in 2SC phase from Ref. [20, [20 and calculated symmetry energy. To estimate the contribution of 
perturbative gluonic interaction, iso-spin dependence of gluon rest masses is calculated by using high density effective 
theory (HDET) [Mil. 


This paper is organized as follows: in Sec. II, brief summary for quark and gluon modification by HDL resummation, 
the thermodynamic potential and the symmetry energy in the normal phase are presented. In Sec. Ill, the gluon 
rest masses in iso-spin asymmetric 2SC phase is calculated by following HDET and its contribution to the symmetry 
energy is estimated. Discussion and conclusions are given in Sec. IV. 
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II. SYMMETRY ENERGY AT LOW TEMPERATURE 
A. Nuclear symmetry energy from equation of state 

Discussion for symmetry energy can start from a finite nuclei with A nucleons Q. The Bethe-Weizsaker formula 
for the nuclear binding energy is given as 


mtot = Nm n + Zrrip - Eg/c 2 , 

E B = ay A — agA^ — ac(Z(Z — 1))A 3 — aAl 2 A + 8(A, Z ), (1) 

where I = (TV — Z)/A. The fourth term in Eg accounts for the total shifted energy due to the neutron number excess. 
The coefficient of fourth term a a represents the shifted energy per nucleon. 

When the atomic number A becomes large, one can approach a limit where parameters describing the system 
become continuous and a mean field type of approach appropriate. In such a system, the symmetry energy can be 
defined from energy per nucleon number: 

E{ - PN a J) = E(p N , I) = E 0 (p N ) + E sym {p N )I 2 + 0(/ 4 ) + • • ■ , (2) 

1 d 2 - 

E S ym(pN) = ^-qJ^E(pn, I), (3) 

where A is the atomic number, p E the nuclear medium density and I the asymmetric parameter I = (TV — Z)/A —> 
[fin — Pp)j{fin + Pp)■ The neutron and proton densities can be found as p n = \pn{^ + I), p P = \pn{^ — I) respectively. 


B. HDL resumed thermodynamic potential 


In this section, we work in the imaginary time formalism (t —^ —it, 0 < t < /3) following Ref. 
function Zq can be obtained in QCD degree of freedom: 


]. The partition 


Zq = Tr exp 


-P(H-p-N) = ip, A, ri) exp 


- j dr J d 3 xC E (ip,ip,A,r]) 


( 4 ) 


where the Euclidean QCD Lagrangian C E is given as 


C E = \f; v F; v + ^A “) 2 + fj a (d 2 S ab + gfatcdpA^ 

+ ^2 \fs d T$f + + m f )4>f - Pftpjipf - g^fMf 

f 


( 5 ) 


where p a is ghost fields and the subscript / represents the quark flavor. From Eq. (U) and Eq. ([5]), the free propagators 
of fields can be obtained as 


Sf(P) = 


mf-f 


— {iu) n + gf) 2 + fP + m 2 ' 


DtyAQ) = 8 


ab 


= 6 


ab 


' 1 

Q 2 

1 

Q 2 


^ QyQv 




L 

[IV 


Q 2 


p T ) 
aiiu 


£ QnQv 

£ QyQv 
Q 2 Q 2 


( 6 ) 


( 7 ) 


where Sf{P ) is the free quark propagator, D'p {Q) the free gluon propagator and £ the gauge fixing term. The 
longitudinal and transverse projection operators are defined as 


P ij = S ij ~ diQj, p L = p Ii = 0 , 

r)L r _ QnQu pT 

Mu' — „ 2 Mii/- 


( 8 ) 
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In the imaginary formalism, the energy integration at each loop changes into a sum over discrete Matusbara frequencies: 



where T^2 is over ioj n = i(2n)irT and iCj n = i(2n + 1)7 iT + fXf for bosons and fermions, respectively. 


1. HDL resumed propagator 

Before calculating the gluonic interaction to the free energy, we summarize the medium modification of the prop¬ 
agator. First, we consider the gluon case. Suppressing light quark masses and soft external momenta (m///x/ — 
0, K M — Qp ~ Kfj) against the hard scale of internal loop (K^ ~ the gluon polarization tensor in the dense 

matter can be calculated as QM1 

n£(Q) = g 2 S ab j l^jTr [ 1im S f {K) 1v S f {K - Q)] 

= m 2 5 ab J ^ (s ijA 6m + , (10) 

where = (— i, k = k/\k\) is light-like four-vector, Q M = (—w, q) is the Euclidean gluon external momentum and to 2 
is defined as 

m 2 = \g 2 T 2 [c A + i»/) + (11) 

v ' / 

At low temperature, only HDL contribution becomes relevant: m 2 => \g 2 Y f kIk 1 ~ ff 2 ff 2 0- The polarization 
tensor can be decomposed into loop structures or polarizations as follows: 

n“t(Q) = 5 ab (n^(Q)| q . h + n^(Q)| q . a ) ( 12 ) 

= <5 ab {SU l (Q)P^ + SU T (Q)P^) . (13) 

First, n u „(CJ)| q „ h and n #iy (Q)| q _ a represents the quark-hole contribution and the quark-antiquark contribution respec¬ 
tively [35l l36(j : 


n^(Q)| q -h = rn 2 J ^-1 + , 

n^(Q)| q . a = m 2 J ^ (^ 4^4 + k^kj). 

Second, the polarization 5H L (Q) and <5n T (<3) can be obtained as follows [33 ]: 


su l (Q) = -2 m 2 ^Q 1 , 


dn T (Q) = m 2 (y~^j 



(14) 

(15) 


(16) 

(17) 


where Qo(x) = \ In \{x + l)/(x — 1)], Qi{x) = xQ o(x) — 1. Here one can find that the longitudinal rest mass comes 
mainly from the quark-hole contribution. 

For the gluon propagator, if the external gluon momentum is soft (Q — g/x), all diagram where HDL is inserted 
as 1PI self energy scales in the same order as the bare gluon propagator (1 /Q 2 ~ 1 /(ff 2 ff 2 ))- The resumed gluon 
propagator is given as follows [34]] : 


*d£(Q) = s ab 


Q 2 


sn L (Q) P ^ 


Q 2 + 5n T (Q) 


p 




£, QnQv \ 

Q 2 Q 2 ) ’ 


(18) 
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FIG. 1. Diagramatic description for ln-Zn, 0 + In Zn g HDL - Each diagram is presented with its symmetric factor. 


Next we consider the quark case. If the external momenta can be neglected against the hard scale of internal loop, 
the quark self energy can be calculated as 0,113 


£(P) = -g 


2 ah 


drK 

, w 

2 nm fc 

47rp.f 


lfi S F (P - K) lv Df v {K) 


(19) 


where = (—u,p) is the Euclidean quark external momentum, r a is the fundamental representation of SU(3) 
generator normalized as Tr[r a r b ] = ^S ab and m 2 is defined as 



( 20 ) 


where Cp = (N 2 — 1)/(27V C ). Again, only HDL contribution becomes relevant at low temperature: m 2 => 
\g 2 C F i±j/ p 2 ~ g 2 g} SI 1 ]- By similar argument with the gluon case, the resumed quark propagator can be given as 
follows [341 ] : 


where 


*S{P) 


1 

f + S(P)’ 


f + E(P) 

A 4 

A s 


A474 + A s 7 • p, 

( 

i \iuj — 7-7-Q0 


rrif 

|p1 + M 


\ 0 \ 

2 r 


IpI 


i-%Qo 

\p\ 



( 21 ) 


( 22 ) 


2. Calculating the thermodynamic potential with HDL resummation 
The free energy can be written as 


= = (23) 

where In corresponds to the sum of all connected ring diagrams. If the interactions are weak enough, one can 
calculate elements in ln^o by perturbative method. However, as the soft excitations in the hot/dense matter have 
non-perturbative loop corrections, this quantity should be obtained by accounting for all equivalent contribution 
through resummation. In the low temperature and dense limit (T ~ gg g), gluonic excitation will be mainly soft 
(Q ~ T ~ gg). In this limit, the physically relevant leading strong interaction contribution comes from the leading 
quark loop (Fig. [[fa)) and from HDL resumed gluon ring diagrams (Fig. |TJb)-(d)). 

The simplest diagram is Fig. Hlb): a soft gluon ring with hard quark loop correction which can also be seen as 
a hard quark ring with soft gluon loop correction. As the gluon line is soft, equivalent HDL corrections should be 
accounted for as a sum of infinite series (Fig. [Hc)-(d)): HDL resummation of gluon ring diagrams. 
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Meanwhile, this type of resummation may not be suitable for quarks as the phy sically relevant quark excitations 
at low temperature and dense limit (\p\ ~ gj) can not have HDL corrections |35[. In this limit, gluon loop should 
carry soft momentum so that the correction will be suppressed by higher orders in a s . However, as there can exist 
suitable intermediate state for HDL resummation among the loop integration of the quark ring trace, we will adopt 
the resummation result reported in Refs. [26lf27j. 

In this work, we consider 2-flavor light quark matter for simplicity. Moreover near the beginning of the quark phase, 
the baryon density ranges from 3-5 po, and the main constituent will be the up and down quarks only. In these region, 
strangeness is not the main excitation mode. Moreover, HDL correction is not relevant for strange quark as its mass 
(m s ~ 95 MeV (H}) is comparable to fj,f. 

The ideal quark gas contribution (Fig. [Ha)) can be calculated as 


]nZn 0 ~/3V 


■ E 

f=u,d 


(24) 


where the absolute divergence known as “vacuum energy” has been neglected. Hereafter, the matter independent 
divergences will be neglected as we are only interested in the finite difference between thermodynamic potentials. 

The ideal gluon gas contribution can be neglected (In-2L 0 — (T ~ g/x) 4 ). The remaining HDL approximated 
gluonic ring diagrams should be resumed. We denote the resummation as “the gluon resummation” (ln.Zo HDL ), 
which corresponds to the diagrams Fig. [Tfb)-(d): 


InZo 


= C 




(N? - 1) 
2 

2 r, 


pv 


pv 


d 4 Q 


(2tt) 

d 4 Q 
(2 tt) 4 


4 ln[l+n M „(Q)D/'*(Q)] 


In 


l + <m L (Q)ij 


2 In 


1 + SU t (Q) — 


where C and T are defined as follows: 

d 4 Q 




t=- ^ n \ pv 


(27T) 4 

d 4 Q 
(2 tt) 4 


In 


In 



„ ? t VjJ , VjJ 

1 + m 2 1 -In — 


1 + 


Next, we take T —> 0 limit and use the spatial d-dimension regularization, 


d 4 Q 

(2n) 4 


du> 

(2?r) 


3-d 


d d q 
(2^)3 ’ 


(25) 

(26) 

(27) 

(28) 


where the discrete sum over Matsubara frequencies T^2 n changes into a finite continuous integration, f dui/(2 tt), in 
the T — >■ 0 limit. One can then regularize the divergence coming from the remaining d = D — 1 spatial dimension with 
g 2 = <?Q/I 2e = <?q/ 2 3_<1 . By extracting out the overall mass dimension as p 3 ~ d in front of the d-dimension momentum 
integration, one can set the coupling constant to be dimensionless and define the one-loop counterterm in lniJQ g HDL . 
A similar rigorous regularization scheme in Coulomb gauge has been already reported in Ref. [25j |. In the present 
work, the calculation is performed in the covariant gauge and the result is totally equivalent to the calculation of 


Ref. [25|, [27( at T = 0. 

As both 5H l (Q) and 5H L (Q) depend only on ioj/q, the integration variable oj can be re-scaled as ui —> qu> to scale 
out the dimensions of the integration: 


£ = - 3 - 2 


duj / dqq d In q 2 + C(cj) 

Jo L 


2 (27r)(27r) d ' 

where the absolute divergence proportional to / 0 °° dqq d In q 2 has been neglected and C(lu) is given as 

W--*_2 A *^ + l 

JL(u>) = m 1 ——In —- 

V 2 ILO — 1 


(29) 


( 30 ) 
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Using the following integration formula [25J 

roo p t l+op p ! l—a \ 

/ dkk a ln(£; 2 + m 2 ) =- - - V 2 / m a+i _(_ absolute divergence, 

Jo a + 1 

and the MS subtraction scheme /z 2 = (47r/e 7 )/2 2 , the integration can be arranged as 

1 dU 3 (to 2 ) 2 


(31) 


£ = (IV 2 - 1)/3U 


PV 


a 


(27r) (27r) 3 4 

2 r 


i i m2 ^ 

1 — in — 7 a — 

7T/U / 


2 \ 1 
a — fd H—a 


V Ui 

Ur Z-^ tt 2 


\f=u,d 

where N c = 3 and the constants are calculated as 


l-t.2 —ln( - 1 ” 


a s \ a- 


finite 


l(u>) = 1 + u (arctanw — — ^ , 


dw/(w) 2 = -tt( 1 - In2) = 0.321336, 


= -0.176945. 


pOO 

P= duil(u) 2 In/(d)) = — 0. 
Jo 


can be obtained in a similar way. After using the same re-scaling u> —> qoj, the integration can be written as 


T=- 


(N 2 - 1 ) i dn d _ 3 _ 


pOO pOO 


-pv 


( 2 ?r) (2ir) d J Q 


/z 3 d 2 duj G?gg d ln q 2 +T(u>) 


where T(w) is given as 


= ± - 


IW , 10J + 1 

_o- 1 -In—- 

ui 2 + 1 2 iuj — 1 


(32) 

(33) 

(34) 

(35) 
s 

(36) 

(37) 


Again, using similar method, the integration can be arranged as 


r = (n 2 - i)pv 


2^2 


1 dVLj, (to 2 ) 

(2tt) (2tt) 3— 8" 


to 2 \ 1 _ 1 5 1 _ 

1 — hr 9 -cr — -p + —a 
2ttm 1 / 2 2 2e 




a, 


’ 7T 1 tj-2 


\f=u,d 


1_ln ( £ _llias ) 5“-^ 


,f=u,d 


finite 


where N c = 3 and the constants are calculated as 

.-,2 


/-\ U / TT _\ 

t(w) = — _ +w - - arctan ui , 
or + 1 V 2 / 

r°° I 

a= dujt(Co) 2 = —7r(—5 + 8 In 2) = 0.142727, 
Jo I 2 

poo 

P= duit(uj) 2 In t{ui) = — 0.2 

Jo 


= -0.200869. 

The overall divergence has been subtracted by the following counterterm (e —► 0) 


AlnZn =-(3 2 - l)/3y 


1 dklo (to 2 ) 2 

(27t) (27t) 3 4 


a+-a 


(38) 

(39) 

(40) 

(41) 

(42) 


For the quarks, the resummation can be done in a similar way as in the gluon case [271 ]. We adopt the result from 
Ref. [13] and denote it as “the quark resununation” (In Zn„ HDL ), which corresponds to the diagrams in Fig. [2] 


In -27) r;iII DL £ ^2 

f=u,d 


-4 


(v) 


8 4 , 

3 - 9 ^ 


(t) : 


(43) 




























+ ••• 



FIG. 2. Diagramatic description for In Zn q HDL - Each diagram is presented with its symmetric factor. 


Finally, lni?n = ln-ZA 0 + In Zq hdl + \i\Zq HDL can be written as an expansion in ct™(lna s ) m as follows: 




f—u,d 


8 4 , 

-7r" 

3 9 


(t)‘ 


2 r 




l ,f=u,a 


1 - Ina* - In [ ^4] )A 1 -A 2 -aln2 

, f=u,d 


(44) 


where the constants Ai = a + 4 a and A 2 = /3 + ^/3. The thermodynamic quantities can be obtained from ln-Z^: 


»(M) {H)-jl-{N) 1 , 

y y ^y ln n 




f-u,a 


,8 4 . 

+ l 3 -r' 


(v) ! 


-5;l E^ 

\f=u,a 

( , (Ni) 13 

M =— = WW* n ’ 


1 —lnct s —lnl Y'' ^-7 —t I 1 Ai — A 2 — a In 2 


1-4 


(tMFH(t) 


E^ls 

\f=u,d 


^ - In a s - In I ^ I I Ai - A 2 - a In 2 


vf=u,d 


M2 


, , (7?) 1/3 1 _ 3 , , 

eW =— = {ap-T'et) “ 


y 

3 


= - y 4 i — 4 f—) 

4 ^ ir 2 V 7 r / 


f=u,d 


8 4 , 

-7b 

3 9 


(t) ; 


+<l E^ 


\f=u,a 


^ — In a s — In I ^ I )Ai-A 2 -aln2 


\f=u,a 


Hi 


where a s (ha) is assigned for 2-loop renomalization in MS scheme: 


(M4) = 


47T 


/3 0 ln (mI/A^) 

Po = (11 N c - 2JV»/3, 

Pi = (347V C 2 - 131V c iV/ + 3Nf/N c )/6, 


1 2ft In (in (^/Af^)) 


(45) 


(46) 


(47) 


(48) 


where h 4 will be parameterized as 1.5/i < m < 4/z. /r is the quark chemical potential at iso-spin symmetric condition. 
As discussed in Ref. [2(| |13, S], if higher order calculation can be obtained, a s can be determined by solving the 
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FIG. 3. (Color online) The ratio P(/x)/Pid ea i(/x) with /i4 = 2p for red dot-dashed and black solid line, and /m = 4 fi for red 
dotted and black dashed line. The upper two lines are obtained with gluon resummation results and the lower two lines from 
adding the quark resummation to the gluon resummation. 


self-consistent gap equation for the thermodynamic potential. On the other hand, since we have calculated only to 
leading order in the HDL perturbation, the effect of a s (^ 4 ) is investigated by varying 1.5/x < /x 4 < 4/x. 

Fig® shows the ratio P(/x)/Pideai(/x) with various /X 4 . The upper two lines are obtained with contributions coming 
from gluon HDL resummation results and the lower two lines from adding the quark resummation to the gluon HDL 
resummation. Red dotdashed and black solid line are obtained with /i 4 = 2/i. Red dotted and black dashed line are 
obtained with /x 4 = 4/x. One can find that HDL resumed P(^x)/Pideai(/x) varies in a large band in Fig® As can be 
expected, P(/x)/Pideai(p) approaches to 1 as /x 4 becomes larger. Here one notes that the gluon resummation gives 
the enhanced ratio (P(/x)/Pideai(/x) ~ 1.1 at yu = 1 GeV), while inclusion of the quark resummation gives the reduced 
ratio (P(/z)/Pideai(M) ~ 0.8 at /x = 1 GeV). The reduction of pressure due to the inclusion of quark contributions 
is consistent with previous studies that include a massive strangeness and the quark interaction part up to 0 (a 2 ) 
HH |H]. Based on the result of Ref. [ 52 J, the authors in Ref. [53j also find a similar behavior of the pQCD ratio 
(0.6 < P(/x) /Pideai(/x) < 0.8 near /x = 1 GeV) as shown in Fig. 2 of their work. While the quark resummation 
enhances the 0 (a s ) correction compared to the weak coupling expansion [26l \ 2 T\, introducing the quark phase in the 
core of a neutron star in general allows for a more massive neutron star that is consistent with recent measurement 
as was discussed in Ref. j53[. To fix the scale, we plotted the quark energy density in Fig. ® One notes that the 
value becomes negative at lower fj, when ^4 = 1.5 fi. This behavior is due to the increase in a s at small scale, which 
is a consequence of asymptotic freedom. To make sure the quark energy density be positive, we choose fj ,4 = 2.5 /x 
throughout this section. 

From Eqs. 65D -6 3 , one can find that the thermodynamic quantities depend mainly on quark and gluon rest mass 
ro 2 ,rn 2 ~ g 2 /x 2 / 7 r 2 . In the HDL resummation, the gluonic contribution is dominated by the longitudinal gluon rest 
mass. Iso-spin asymmetric properties can be accounted for by assigning asymmetric quark chemical potential to the 
quark and gluon rest mass. 


C. Cold dense matter symmetry energy at normal phase 


In the iso-spin asymmetric matter, the light quark chemical potential can be defined as follows: 


Pd=P 




( 49 ) 
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/i( GeV) 

FIG. 4. (Color online) Quark energy density e(/x) with various /m. Blue dot-dashed line is for ideal quark gas. Black lines 
represent the result for quark and gluon resummation. Black dotted line represents the result with = Ap, black dashed line 
with /m = 2.5 fi and black solid line with /m = 1.5 fi, respectively. 



0.36 0.38 0.40 0.42 0.44 0.46 0.48 0.50 

/i(GeV) 


FIG. 5. (Color online) Quark matter symmetry energy at normal phase (/m = 2.5 fi). Blue dot-dashed line represents symmetry 
energy of ideal quark gas. Red dotted line represents symmetry energy with gluon resummation only while black dashed line 
represents the result after adding the quark resummation. 


By assigning this relation to Eqs. gj - grp , quark matter symmetry energy in the normal phase can be obtained from 
energy per baryon number as follows: 

—= E(fi, Ib) = E 0 (fx, Ib) + E sym (p)I^ + 0(1 %) H- , 

Pb(pOb) 

1 d 2 - 

E S ym (/t) = 2 ^ gj2 Ib )j 

where p B = {p u + Pd)/ 3, Ib = 3 (p u - Pd)/{p u + Pd)- 

The quark matter symmetry energy is plotted in Fig. [5] Here one can check that the symmetry energy is reduced 
by HDL resummation compared to the ideal quark gas case. It reduces slightly when the gluon HDL resummation is 
included and more when the quark resummation is added. Hence, one can conclude that HDL resumed interaction 
makes it easier for the cold quark matter to be iso-spin asymmetric: it costs less energy to reach the same iso¬ 
spin asymmetry in comparison with the ideal quark gas. When temperature becomes even lower, additional non- 
perturbative effect should be considered. We discuss an effect in the next section. 


(50) 

(51) 
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FIG. 6. (Color online) Ratio P(g) / P\dea.\{p) with various A. Black dashed line represents P{p) / Pidea,i(g) with A = 150 MeV. 
Black dotted line represents P{p,)/Pideai(g) with A = 50 MeV. Red dot-dashed line represents P (g) / Pidea\(p) of the gluon 
resummation with fM = 2.5 g at normal phase as a criterion. 


III. SYMMETRY ENERGY AT EXTREMELY LOW TEMPERATURE 


Now we consider the extremely cold matter case (T -C gg). As we have found in the previous section, the gluonic 
contribution to symmetry energy depends mainly on the iso-spin asymmetric quark and gluon rest masses which 
comes from HDL resummation. In this section we will concentrate on the iso-spin asymmetric corrections which come 
from extremely cold matter. 

In such a cold condition, for ideal quark gas, all fermion will be confined in their Fermi sea without large fluctuation. 
When interaction turns on, one can consider nontrivial correlation between fermions and related effect. According 
to effective descriptions based on Wilson’s renormalization group approach 0NI3 , when one scales the momentum 
to near the Fermi surface, the interaction between fermions with opposite momenta becomes marginal. So if the 
interaction is attractive, it is natural to form a bosonic condensate for two fermions on the same Fermi sea with 
opposite Fermi velocity (BCS pairing) [18]. In QCD, one-gluon exchange via color antisymmetric channel naturally 
provides attractive interaction. So in the extremely low temperature limit, one should consider the effects of color 
superconductivity 01. 


A. 2-color superconductivity (2SC) and symmetry energy 


First, we confine the density region to be similar to that of previous section: 3-5po- The corresponding quark 
chemical potential then ranges from 0.35 GeV < g < 0.5 GeV. In this region, there is a mismatch of Fermi momentum 
between light quark and strange quark (p F ~ p F ~ g p s F = gl — fn 2 s ) even if same numbers of light quarks and 
strange quarks exist. Thus when str ang eness is not the main excitation, only light quark flavor can participate in 
formation of BCS condensate (2SC) [24j. So again, we consider only light quark flavors in this section. The spin-0 
channel forms the most stable diquark condensate as quark pairs from the whole iso-tropic Fermi surface can contribute 
to the pairing. This type of condensation requires the quark pair to be on the same Fermi surface but with opposite 
velocity and spin alignment, such that they are in same hclicity state. As color configuration is already antisymmetric, 
to obtain total antisymmetric wave function, quark flavors should be in the antisymmetric configuration. The diquark 
condensate can then be summarized as follows fiEI [33}: 

{i>L, a % C ^L,Pj) = - = yea/33Gj3, (52) 

where C = icr 2 , A is the superconducting gap size, the Greek index represents color, the Latin index represents 
quark flavor and L , R represents chirality of the quarks. The condensate can be included in the interaction terms of 
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Lagrangian as the invariant coupling [32| [33j : 






(L —> R) + h.c., 


(53) 


where e = ia 2 on SU(2) color-flavor index. 

Although the existence of the 2SC phase is based on the iso-spin symmetric Fermi sea of the light quarks at zero 
temperature, the phase can persist up to some critical temperature and iso-spin asymmetry determined by the Clogston 
limit (5/i = A/\/2 [2(1 Rot . Specifically, the PNJL model calculation [50] suggests that the critical temperature is in 
the order of 50 MeV. Moreover, even if the pure neutron matter (7g = 1) changes into quark matter, the difference 
between Fermi sea of u and d quark will be Sg ~ 0.23/1. Hence, assuming that the gap also grows with the chemical 
potential, the Clogston limit can be satisfied at higher density. In realistic heavy ion collisions involving ordinary 
dense matter ( Ib < 1) the iso-spin difference will be smaller (Sg < 0.23/i) and and the possibility to probe 2SC matter 
remains. Such a matter can be produced in the planed FAIR experiment [5l| . 

If one considers only the ideal quark contribution and the correction coming from the non-perturbative gap contri¬ 
bution, the free energy can be written as follows [H, |4(j : 


^a(p) — ^g,o(p) A ; a(a0 


1 Nc 


,,4 2SC 2 


f=u,d 




N c 


_y 

_ 12 2—^ ^2 


f=u,d. 


fff 

7T 2 


2 A 2 


47T 2 


(54) 


where the approximation in the last line follows from Ref. [2(1 [29[. In the 2SC phase, the 2-color and 2-flavor states 
participate in the condensation and hence share the same chemical potential. Here, we ignore u and d quark chemical 
potential asymmetry coming from the difference in their charge. Then, the correction from 2SC can be written as 
(5H 9 .a(m) = —4(/x 2 A 2 /47t 2 ) without the iso-spin index. 


The gap size A has been determined in several studies. The first work in non-perturbative regime was reported bv D. 
T. Son as A / g = (b/g 5 )ex p (~3tt 2 /y/2g) by solving the self-consistent gap equation in using HDL approximation 211. 
Schafer and Wilczek [23 determined the constant to be b = 2507r 4 and the corresponding gap size in g = 400 MeV 
to be in the order of 100 MeV. The estimation from HDET [33j was found to be in the order of 50 MeV, while model 
calculation [23J gives the value in the order of 100 MeV. 


In Fig. [6] we plotted P(g )/ -Pideai(p) in the range of 50 MeV < A < 200 MeV. As one can expect, as A becomes 
larger, the pressure deviates more from the ideal quark gas limit. When A = 150 MeV, P(g )/ -Pideai(p) shows similar 
behavior to P(g) /Pideai(p) at normal phase (with only the gluon resummation, g^ = 2.5/x). Also with the gap size 
A = 150 MeV, large enough Clogston limit (Sg ~ 0.23 g < A/-\/2) can be obtained in 3-5po region. For gluonic terms, 
at 2SC phase, HDL effective Lagrangian can not be inserted directly to the present Lagrangian as the quasi-loop 
correction in 2SC phase is different with the HDL of normal phase. However, it was argued in Ref. J4l| that the 
coupling in such cases will be weak, hence, may give minimal effect to the free energy. In this section, we choose the 
A = 150 MeV to match the value of the free energy calculated in both phases as shown in the upper two lines in 
Fig. [6] 


The symmetry energy in 2SC phase can be obtained in similar way from the energy per quasi-baryon number: 


= E fa T b) = E o(P Ib) + EsymWl + 0(4) + • • • , 

Esymif 1 ) = qj2 ^ b)j 


(55) 

(56) 


where p/j = (p U npaired + Ppaired)/3 and /g = /g/3 as asymmetrizable unpaired states reduce to 1/3 of the normal 
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FIG. 7. (Color online) Quark matter symmetry energy at 2SC phase. Blue dot-dashed line represents symmetry energy of the 
ideal quark gas. Black solid line represents symmetry energy at 2SC phase with A = 200 MeV. Black dashed line represents 
symmetry energy at 2SC phase with A = 150 MeV. 



FIG. 8. Relevant one loop diagrams for the gluon rest masses, (a) represents the ordinary polarization part, 
the momentum independent tadpole part which gives counterterm to transverse gluon mass. Double line in 
propagation in Dirac sea mode I30l. l3lll . 


(b) represents 
(b) represents 


phase. From the quasi-quark number density and the energy density can be obtained as follows: 


Pi((J>) 

PAi(p) 

ea(m) 


P, 


A 2 


1 M? 

3 7r 2 

iMf _ 

3 7T 2 ^ 2-7T 2 

^unpaired (aO ' 
unpaired a 

i ^ ft 

4 J 7r 2 


^paired (m) 

, 2SC r 4 
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(57) 

(58) 


(59) 


where fi% = [i( 1 =F Ig) 3 is assigned only for the unpaired states, while the other quasi-quark states share symmetric 
chemical potential /r. The symmetry energy in 2SC phase is plotted in Fig. [7] It is substantially larger than the 
symmetry energy in the normal phase: it costs more energy to reach the same iso-spin asymmetry in comparison 
with the normal phase. This large difference may lead to quite different phenomenological prediction, which will be 
discussed in last section. As the gluonic coupling will be weak [41], the result in Fig.[7]could be dominant contribution. 


B. Gluon rest masses in iso-spin asymmetric 2SC matter 

While we expect the gluonic contribution to the symmetry energy in the 2SC phase to be weak, it is worth while to 
estimate it. As we have found in the previous section, contribution of HDL resummation to the symmetry energy in the 
normal phase is determined dominantly by the quark and gluon rest masses. Among them, the gluon resummation 
with hard quark loop correction seems suitable for physically relevant condition. In the 2SC phase, the gluonic 
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a,b = 1, 2, 3 

a, b = 4, 5, 6, 7 

a, b = 8 

Paired (A=0) 

s-ab 1 cab 

SAA = 4 ° 

tab 1 cab 

SAA 8^ 

cab 1 cab 

SAA 12 ^ 

Paired (A=l,2,3) 

= \s ab 

Ga = 

Ga = P 2 5 ab 

Unpaired (A=4,5) 

Ga = 0 

Ga = 

Ga = 


TABLE I. List of Ga- Small letters a,b represent adjoint color index of the external gluon. Capital letter A represents the 
quasi-quark state in the loop of the tadpole diagram (Fig. [Sfb)). 


excitation will be soft or super-soft (Q ~ T <C gg) and we assume that the qualitative behavior of the symmetry 
energy in 2SC phase also depends mainly on the quasi-quark loop correction in small momentum limit. By using 
HDET 3(]- 33|, the rest masses of gluon can be simply calculated. In HDET, matter part of Lagrangian can be written 
as follows 3(1 Ol: 


= + Ai/7o)VT 

/ 

= E E • D ^f+ - ^/+ 2m/ + iv . D {P±)2 ^ f+ ) ’ (60) 


where = (1, Vf) and V M = (1, —v/) with Fermi velocity Vf , p± = 7 with 7^ = 7^ — 7^ and 71 1 = (7 0 , VfVf 7) 
and i/jf+ represents positive energy projection of quark field. The irrelevant quark mode ipf- is integrated out by 
using equation of motion: ipf- = —*7°/ (2/// + iV ■ D)p±ipf + . 

In 2SC phase, the quasi-quark states can be represented as follows [32| : 


/ (Aa)cm iA A 

i>± ,m = 2^ - 7k ^±> X = 


lf>+ 


A=0 


V2 


Cip A 


A * 


~ 1 2 ~ 1 1 

Aq = -^=As + —I', A a = Aa {A = 1, 2 ,3); A 4 = —^=(A 4 — A 5 = —j={\q — i\ 7 ), 


(61) 


where A = 0,1, 2,3 denotes gapped states, A = 4, 5 denotes ungapped states and A,; is Gcll-mann matrix. Here, ‘±’ 
do not represent the energy eigenstate. It represents direction of Fermi velocity. Incorporating this representation 
and the gap interaction term (1531) . the Lagrangian can be written in Nambu-Gorkov form [32, ,33| as 


C = 


_ -F a F a » v + V V l^t ( lV ' 95 An Aab \ Y B +iaA a ( iVI " KAaB 0 ^ 

4 ' VxfeoL V ^ iV-dS AB ) X + ' SA “ X ( 0 ->VV A , B j 

/ 1 P cd D 

+ g 2 AlAi x ^ I *rW.D*AB 


X 




2 p f -\-iV • D' 


+ (L — y i?), 


(62) 


where P^ = g^ v — \(V ll V 1 ' + V V V^) and the constants are defined as 

KAaB = ^Tr[\ A T a \B], £aB = iTr[AAT c T d A B ], 

A A b = y Tr[ea A ea B ] (A , B = 0 , 1 , 2 ,3, otherwise A A b = 0 ), 

= A a Sab with A a = (-A, A, A, A, 0,0), (63) 


where r a is the fundamental representation of SU(3) generator. 

For gluon self energy, the diagrams in Fig. [ 8 ] are relevant. The diagram Fig. [U(a) can be calculated from ordinary 
current-current correlator 


n 




half 

d A xe~ iqx J2 

Vf 


(n\T[j;(x,v f )Jt(o,v f )]\n), 


(64) 
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a,b = 1,2,2, 

a, 6 = 4, 5, 6, 7 

a, b = 8 

ngo(o) 

0 

1 „_2 cab 

2 m 0 

m 2 8 ab 

— n “/ (o) 

0 

ik 2 £/V/A 2 )M ab 

\[-m 2 +g 2 Y,T^)h r 2 )]M“ 6 


TABLE II. Debye and Meissner masses in the asymmetric 2SC matter. Small letters a,b represent adjoint color index of the 
external gluon, m 2 = (g 2 g , 2 /tv 2 ). 
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FIG. 9. (Color online) Symmetry energy as a function of the Fermi momentum. Left side: symmetry energy in the hadron 
phase calculated by using QCD sum rule [46}. Blue dot-dashed line represents the kinetic part of symmetry energy of the ideal 
nucleon gas. Red dotted line represents the kinetic part of interacting matter. Black dashed line represents the total symmetry 
energy (kinetic + potential). Right side: quark matter symmetry energy in the normal and 2SC phase. Blue dot-dashed line 
represents symmetry energy of ideal quark gas. Red dotted line represents symmetry energy in the normal phase with HDL 
resummation. Black dashed line represents symmetry energy in the 2SC phase. The phase boundary has been set artificially. 


where |f2) denotes 2SC ground state and J£{x, Vf) can be obtained from the second term in the matter Lagrangian. The 
diagram Fig. [5]}b) directly comes from the last term in the matter Lagrangian which has explicit iso-spin dependence. 
In iso-spin symmetric condition, with counterterm dll“ h (0) = — i m 2 5ijS ab which should come from Fig. (5]b) type 
diagrams as it is transverse component, the gluon rest masses can be matched with Debye and Meissner masses 
calculated in Ref. [36}. In iso-spin asymmetric condition, only thorough this counterterm, iso-spin asymmetric part 
of the rest masses can be accounted. 

The iso-spin asymmetric part in the whole counterterm 5II“ & (0) can be determined by calculating Caa- With 
in TableU the Debye and Meissner masses in the asymmetric 2SC matter can be obtained. The Debye and Meissner 
masses are listed in Table dn From Tabic HU one can find that the rest masses of the gluon in adjoint color a = 1, 2,3 
vanishes. As this type of gluon carries the unbroken SU(2) color charges which participate in the 2SC condensation, 
the rest masses vanish trivially. The gluon in adjoint color a = 4, 5,6, 7,8 has rest masses. It can be noted that 
only the Meissner mass is iso-spin weakly dependent. For the static quantity such as thermodynamic potential, the 
contribution of the Debye mass is most important. Consequently, in the 2SC phase, the contribution of the Meissner 
mass to the symmetry energy will be minimal. However, for dynamical quantities, the situation can be different. 
The spatial transverse mode of super-soft gluon can acquire nontrivial loop corrections that can be as large as HDL 
|42l - l45| . This means that dynamical process which accompanies super-soft gluon in 2SC phase may strongly depend 
on iso-spin asymmetry. Thus such a process may critically influence observables that can discriminate the existence 
of 2SC phase and identify the symmetry energy properties of the cold dense matter. 


IV. DISCUSSION AND CONCLUSIONS 


In this work, we calculated the thermodynamic potential and the symmetry energy at the cold dense limit. For the 
normal phase at low but non-negligible temperature, we first calculated the free energy by using HDL resummation, 
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which corresponds to calculating quark and gluon ring diagrams with HDL self energy. From a theoretical point of 
view, the gluon resummation seems to be more suitable for the physical low energy excitations in the cold dense 
matter. For the quark resummation, we optionally considered the effects of the resummation. All the thermodynamic 
properties, calculated from the thermodynamic potential, strongly depends on the quark and gluon rest mass. In 
the gluon self energy as well as the free energy itself, the largest contribution comes from the Debye mass of the 
longitudinal mode, which comes mainly from the quark-hole loop in HDL self energy, whereas the contribution from 
the transverse mode, which comes from quark anti-quark loop, is comparatively unimportant. The symmetry energy 
also comes mainly from the iso-spin dependence of the Debye mass. HDL resummed gluonic interaction gives reducing 
effect to the symmetry energy. 

For the 2SC phase, which can appear when temperature becomes extremely l ow, the free energy containing non- 
perturbative condensation should be considered. Following previous studies [28l, 129}, the free energy can be written 
in a simplified form as given in Eq. man, from which one can extract the symmetry energy in the 2SC phase by 
reducing the formula to have the iso-spin asymmetric factor IB ' The symmetry energy becomes almost 3 times that 
in the normal phase, because 2-color and 2-flavor states are locked in their common Fermi sea, reducing the number 
of available quarks that can contribute to the asymmetry. Although the symmetry energy becomes larger than the 
typical value in the normal quark phase, it still remains smaller compared to the expected value in the hadronic phase 
extrapolated to the phase boundary (Fig- HJ- We also estimated the gluonic contribution to the symmetry energy in 
the 2SC phase. By using HDET j.'iol |.3.lf . we found that only the Meissner mass in the 2SC phase has a small iso-spin 
dependence (Table. [II]). In another words, as far as the gluons are concerned, the Debye mass is important for the 
static quantities, and Meissner mass contributes to the symmetry energy, but with an even smaller reducing effect 
than in the normal phase. Hence, one can just consider the effect of BCS condensation only. 

With the results on the cold dense matter symmetry energy in hand, one can imagine a phenomenological scenario 
in which one can create a cold dense matter through a heavy ion collision. In Fig. [9] one finds that the nuclear 
symmetry energy in the hadronic phase rises stiffly. The nuclear symmetry energy was calculated by using QCD 
sum rules, based mainly on linear density approximation [46| . The OPE terms and formula for plotting density the 
behavior can be found in Ref. (46| . Similar stiff density behavior of the symmetry energy was reported before in 
Ref. including iso-spin dependent interaction channel (NLpd). If the symmetry energy increases stiffly at high 
density region, the neutron sea will become unstable. Then, it may lead to nn —l pA~ type scattering. As the 
unstable neutrons are scattered out, the iso-spin asymmetry of remaining nuclear matter will be reduced. But, as 
proposed in Ref. M, if the mixed phase with quark and hadron appears, the iso-spin asymmetry can remain high 
due to the small symmetry energy of the quark phase: iso-spin distillation, by forming deconfined quark drop in the 
neutron. This distillation effect will cause enhanced n~ yields compared to the case where there is only hadronic 
phase. Our calculation shows enhanced quark matter symmetry energy when 2SC phase appears. So the distillation 
effect will be reduced and the ratio 7r - /7r + will be smaller than the case where only the normal phase exists [5s, 6|. 
Our consideration is in agreement with Ref. Q. If the heavy ion collision experiments can be arranged to produce 
different temperature of dense matter, these two trend can be identified: one goes to the normal quark matter with 
low temperature and the other to the 2SC quark matter with extremely low temperature. The two experiment may 
give enhanced ratio compared to the hadronic phase, but with a quite different value. Such experiment may 

be available at FAIR [5li ] . 

Although we have neglected strangeness degree of freedom for simplicity, considering its effects may give us a 
fruitful insight. Even if they are not the main excitation in the matter, one can expect high multiplicity of the hadron 
excitation containing strangeness. Measuring kaon multiplicity may be good probe for the cold dense matter as it 
has higher threshold energy and do not strongly interact with surrounding nuclear matter (4?l l48j . Also the iso-spin 
dependence of hyperon mass and its interaction scattering cross section at high density regime can be an interesting 
issue. Describing these physics using explicit QCD quantum number may provide new understandings for high density 
iso-spin asymmetric matter. Related study is now in progress. 
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